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搞 要 : 本 文 利用 Kronecker 积 的 性 质 ， 得 到 了 非 线性 矩阵 方程 X = Q 一 A" (Im 8 X — C) -4 存在 正定 

解 的 充分 必要 条 件 。 运 用 有 界 序列 的 收敛 原理 ， 给 出 了 求解 方程 的 不 动 点 迭代 与 无 逆 和 迭代 两 种 移 

代 方 法 。 数 值 例子 验证 了 这 两 种 迭代 方法 是 行 之 有 效 的 。 
XS. 非 线性 矩阵 方程 ， 正 定 解 ， 不 动 点 迭代 ; EER 


分 类 号 : AMS(2000) 15A24; 65H05 中 图 分 类 号 : O241.7 文献 标识 码 : A 
1 引言 
本 文 研究 非 线 性 矩阵 方程 
X =Q- A (oe XC) A (1) 


在 矩阵 集合 Xe p(n) = (X | Im X > C) 内 的 正定 解 ， 其 中 @ 是 m 阶 正定 矩阵 ，C 是 mn 阶 
PERHE, HG QC > 0，4 是 mn xzm 阶 复 矩 阵 。 形 如 (1) 的 非 线性 矩阵 方程 在 控制 理 
论 、 梯 形 网 络 、 动 态 规划 、 统 计 学 等 领域 中 有 广泛 的 应 用 [9]. 


2 EDE (1) 有 和 解 的 条 件 
方程 (1) 等 价 于 
Im 8 X + (Im Q A*)| Im 8 (Im 8 X -C) (In®A)= In QQ. 
x z z 
Y-Il,69X-C, Q-I,6Q-C»0, A=In®4, 


则 方程 (1) 可 转化 为 
Y + A' (I, 6Y !)A-Q. (2) 

将 4= In A 分 块 为 
A= (AT, Asl, PU Ast). 


? 
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HH A = eiT QA, i=1,2,… ,m， 再 令 
则 方程 (2) 又 可 转化 为 


引 理 2.1 方程 (1) 有 解 等 价 于 方程 (3) 有 解 。 
引 理 2.26] 方程 (3) 有 和 解 的 充分 必要 条 件 是 A; 有 以 下 分 解 


A; = W*Zi, $-1,2,-- M, (4) 


Xp W E ARENA (WT, ZT,..., ZD)T ENER. CERET, Y = WW 是 方 
FE (3) 的 解 。 
定理 2.1 方程 (1) 有 解 的 充分 必要 条 件 是 A4 有 分 解 4 = W*Z, Hop W 是 非 奇 异 矩阵 且 有 


W*W + Im 6 (Z*Z) =Q. (5) 


在 此 条 件 下 ，X =Q- ZZ 是 方程 (1) 的 解 。 
证 明 必要 性 : 设 方程 (3) 有 解 立 。 由 引 理 2.2 可 知 ， 有 分 解 式 (4)， 且 有 


(W*, Zi COE Z;)(WT, Zt SD 27) = Imn, 


4W-WàQs:, WaR mW ERR. B 
W-*(eT 8A) 2 e (W™A) = ZiQ*. (6) 
将 (6) 式 的 右边 丈 G 主 分 成 与 (6) 式 中 间 相 同 的 1x m 块 形式 ， 设 其 中 第 i 块 为 Z;， 则 由 (6) 式 得 
ZQ? =eiT & Zi = eT 8(W™*A), 
HZ, = W-*A, i=1,2 m. A W-"A- Z. 于 是 4= W*Z, Zi = (eit ® Z)Q-3, 


&5W-WQ-iü 


Q-iw*wQ-? 十 5 [Q- e; Q Z*) (e; e z)Q i ] 


i=} 


= qiwwa-i «o eene no? 


i=1 
-Qi [w*w 十 ( > ees) & (2*z)|o7* 
i=l 


= Q3G[WWw + Im ® (2*2)] 0 3 = Imn, 
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即 有 W*W + Im 8 (Z*Z) = ÑR. 
充分 性 : SA-W'"Z. W 非 奇 异 ， 且 条 件 (5) 成立， 则 


A = Q2 (eT e 49-3 =Q [aT e (W*2)] 9-3 = 0 iW (eT e Z). 


4 W-WQ-i, Zi 2 (NT 8 Z)Q-?. WA; = W*Z,, WW 非 奇 异 。 结 合 条件 (5) 不 难 验证 


W*W + Y Zt Zi = Imn, 
{=f 
Bl (WT, Z7,- ,2ZT)7 了 是 列 正 交 的 。 由 引 理 2.1 可 知 ， 方 程 (3) 有 和解 ， 从 而 方程 (1) 有 解 ， 
HY = WW = Q-3W*wQ-3$ 是 方程 (3) 的 一 个 解 。 XY = QÝ? = W*W， 代 入 方 
程 (1) 得 
X =Q- AYA =Q- AWWA =Q- Z*Z. 


3 求解 矩阵 方程 (1) 的 数值 方法 


ATXJfS., FOOUHXGeRIAGX. 
算法 3.1 


Xo = Q, 
Xa4412Q-—A'(X,—C)A, n=0,1,2,.…. 


为 了 讨论 算法 3.1 的 收敛 性 ， 首 先 给 出 如 下 引 理 : 

引 理 3.1 设 X 为 矩阵 方程 (1 在 o(m) = (X| X > C} 内 的 任意 解 ， 则 X < Q。 

313.24 WRA > B > 0， 则 当 a € (0,1, 有 4c > B^ > 0; 当 ae[-10) 时 ， 
有 Ba > A9 > 0. 

定理 3.1 WX € y(n) 为 方程 (1) 的 任意 解 ， 则 由 算法 3.1 得 到 的 序列 {X%} 收 剑 到 方 
程 (1) 在 p(n) 内 的 极 大 解 Xz。 

证 明 采用 归纳 法 可 证 序列 {Xh} 有 下 界 针 且 序 列 {X%} AER, MOX) Sx. X X E 
HEO) Egeln) 内 的 任意 解 ， 因 此 {X} 收敛 到 方程 (1) Æ pln) 内 的 极 大 解 Xr. 

下 面 ， 我 们 给 出 求解 矩阵 方程 (1) 的 无 逆 迭 代 算 法 。 不 失 一 般 性 ， 下 文 我 们 假定 方 
EOF = ， 则 条 件 @ > CERA >C 

算法 3.2 


Xo—1,, Yo—(i,—-C)-!, 
Xk+1 = In — A*Y4A, 
Yet1 = Yel2fn — (Rr — O)Yk]), k-012.-. 
定理 3.2 WX, Ynn = 0,1,2,… 是 由 算法 3.2 产 生 的 迭代 序列 。 若 方程 (1) 有 一 个 


解 X € p(n)， 则 方程 (1) 在 y(n) 内 存在 最 大 解 Xr， 且 {Xn} 单 调 递减 收敛 到 Xr。 
证 明 此 定理 证 明 类 似 于 文献 3] 中 定理 3.3， 故 略 去 。 
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4 数值 例子 
我 们 先 来 考虑 算法 3.1 与 算法 3.2 的 计算 复杂 度 。 不 妨 假设 所 有 的 矩阵 都 是 实 的 。 对 于 算 


法 3.1 中 的 矩阵 求 逆 运 算 ， 一 般 采 用 Cholesky 482, Bl X, — C = LLT7， 然 后 再 求解 三 角形 和 矩 
EDF LB =A, RAWA Xn = Q 一 BTB。 容 易 得 到 算法 3.1 每 迭代 一 步 需 要 





nar + m?n? + mn? 
个 Hops， 计 算 中 用 到 了 矩阵 的 对 称 性 。 类 似 地 ， 可 以 知道 算法 3.2 每 迭代 一 步 需 要 3(mm)3 十 
2m?n? + mn? 个 Hops。 不 难看 出 ， 虽 然 算法 3.2 避 免 了 矩阵 的 求 逆 运 算 ， 但 是 其 每 欠 代 步 计 算 
所 需 fops 比 算法 3.1 要 多 。 
下 面 用 数值 例子 来 验证 算法 3.1 与 算法 3.2。 以 下 所 有 结果 都 是 在 Matlab 7.0 中 运行 得 到 的 。 
例 1 考虑 矩阵 方程 


X-Q-A'(I,9X-C) A, (7) 
其 中 
1 2 1000 
01 32 13 9 0100 
m = 2, A= E T ? 一 
5 6 9 244 0000 
5 9 0000 


若 用 算法 3.1 求解 方程 (7)， 经 24 步 迭代 后 得 到 方程 (7) 的 极 大 解 


10.4901 5.9473 
X x Xay = , 
5.9473 235.9229 


AiR, RÆ 
| RGGA)|| e = || X24 — Q + A* (15, 8 X24 — C)TA]| p = 7.0072 x 10715. 


若 用 算法 3.2 求 解 方 程 (7)， 经 49 步 迭代 后 得 到 方程 (7) 的 极 大 解 


10.4901 5.9473 
Xm X49 = , 
5.9473 235.9229 


经 计算 ， 残 差 
| RGGo) || = |lXao — Q + A* (Im 8 Xa9 — C) Al = 6.7875 x 1075. 


由 例 1 可 知 ， 用 算法 3.1 和 算法 3.2 求 解 矩 阵 方程 (1) 是 可 行 的 ， 而 且 也 是 有 效 的 。 
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Positive Definite Solutions of the Matrix Equation 
X-2Q-A'(I,.8X—C)!A 
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Abstract: The nonlinear matrix equation X = Q — A* (Im & X — C)! A is considered in this paper. 
A sufficient and necessary condition for the existence of positive definite solutions to this equation 
is presented by using properties of the Kronecker product. Two iterative methods for solving the 
equation are constructed by making use of the principle for the bounded sequence. The feasibility and 
effectiveness of the iterative methods are verified by a numerical example. 

Keywords: nonlinear matrix equation; positive definite solution; fixed point iteration; inverse-free 
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